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Abstract 


Prediction of Ae 3 temperature in case of steels is quite important both from the practical 
as well as theoretical standpoint. Although, in the past several attempts have been made to 
develop a model to predict the Ae 3 temperature from the composition of the alloys, not a 
single model is known to be full-proof. Therefore, there still remains the need to establish a 
model that will predict the Ae 3 temperature as a function of alloy composition as accurately 
as possible. To this end, the present work is an attempt to model the Ae 3 temperature of low 
carbon steels based on composition. The alloying elements considered are C, Mn, Si, Ni, Cr, 
Mo, and Cu. In all, 182 binary and multicomponent data on Ae 3 temperature are considered 
from different sources. Seven nonlinear binary equations are developed. The predicted Ae 3 
temperature for the binary systems is quite close to the experimental value (residue < ± 1 ° 
C). In case of the multicomponent system, a polynomial model is assumed. This model is 
tested in the following way : keeping the constant term along with the coefficient of the 
binary terms fixed to that of previously established binary models, the coefficients of the 
multicomponent terms are evaluated by stepwise linear regression analysis. By doing so, 
the total error in predicting the Ae 3 temperature is reduced as compared to other available 
models. The standard error in the present model is 7.34, which is much less compared to 
earlier models. The model predicts the Ae 3 temperature over a wide range of compositions. 



Chapter 1 


INTRODUCTION 


Steel is an important alloy of iron having multipurpose usage. It is used in most of the engi- 
neering applications because of its unique properties like high strength, hardness, toughness, 
ductility, and corrosion resistance. It has wide applicability owing to its availability and the 
cost factor. The properties of a typical steel depend on its composition, micro-structure and 
various heat treatments given to it. Although these processing parameters are interrelated, 
some of the desired properties can be obtained in a steel purely by proper heat treatment. 
For this purpose, it is essential to ascertain a priori the austenite to ferrite transformation 
temperature or the Ae 3 temperature. The Ae 3 temperature varies with compositions and 
alloying elements present in a steel. 

Ae 3 temperature of a steel can be determined experimentally. In industries, steels with 
wide range of compositions and components are used. Accordingly, the Ae 3 temperature 
varies. Determination of Ae3 temperature of each and every steel by experimental means is 
very inefficient both in terms of cost and time. Many a time it is not possible to carry out 
experiments. Therefore, efforts have been made to find out an empirical relationship between 
Ae 3 temperature and compositions of steel. To this end, there are several mathematical 
models available in the literature but none of them give satisfactory results. The errors in 
these models are either due to wrong assumptions inherent in the models or the data used 
to develop the model were erroneous. In some models, the alloying elements considered 
were very few while in others, the interaction between various alloying elements were not 
considered. Although the most recent model based on genetic algorithm predicts the Ae 3 
temperature satisfactorily, it still suffers from its limited applicability. This is primarily due 
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to the use of binary models that were calculated over a narrow range of compositions. 

The present work was carried out to develop a mathematical model to predict the Ae% 
temperature in steels based on compositions. In this model, unlike in the previous models, 
a large number of data were considered. In particular, the number of binary data used in 
this model is larger than the previous models. Using these binary data, the effect of seven 
alloying elements namely C, Mn, Si, Ni, Cr, Mo, and Cu on Ae 3 temperature was studied. 
Respective binary iron-alloys models were developed over a wide range of compositions and 
all of them were found to be nonlinear. These binary models were adapted in a polynomial 
model for the multicomponent system and the important multicomponent terms and their 
coefficients were determined by stepwise linear regression analysis. The results obtained in 
the present model were statistically analyzed and found to be more accurate compared to 
existing models in the literature. 



Chapter 2 


LITERATURE REVIEW 


2.1 Binary Phase Diagrams of Iron 

A majority of low carbon steels have two- phase structures, namely ferrite and cementite, 
at room temperature. Heat treatments given to a steel determine the micro-structures and 
phases present in it, at room temperature. Many properties, like yield strength, ductility, 
and toughness, are dependent on the micro-structures and phases present in the steel. The 
knowledge of Ae 3 and Ae 1 temperature of a steel is very helpful for heat treatment. The 
. 4 e 3 temperature represents the temperature above which a hypoeutectoid steel shows fully 
austenitic structure. Similarly, Ae 1 temperature is that temperature where austenite begins 
to form in a hypoeutectoid steel. 

Data on Ae 3 temperature were carefully selected, as in most of the cases the tempera- 
ture refers to either Ac 3 or Ar 3 temperature. Ac 3 temperature represents the temperature 
above which fully austenitic structure is obtained while heating and Ar 3 temperature is the 
temperature below which ferrite begins to form while cooling a hypoeutectoid steel. Thus, 
Ac 3 temperature lies above Ar 3 temperature, and Ae 3 temperature lies in between these two 
temperatures. However, it is not necessarily the arithmetic mean of the two [1]. 

In steels, Fe combines with elements like C, Mo, Si, Cr, and Ni to form commercially 
important alloys. The alloying elements are classified into two types depending on their 
effect on phase fields. They are - 

1. a stabilizers - These alloying elements expand the a-phase field, e.g., Cr, Mo, and Si. 

2. 7 stabilizers - The solutes which expand the 7-phase field are known as 7 stabilizers, 


3 
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e.g., Ni, Mn, C, and Co. 

Depending on the influence of alloying elements, the phase diagrams are divided into four 
categories : 

• Type 1 Open 7 field - The alloying elements like Mn, Ni, and Co, expand austenite 
region towards room temperature and raise the two phase (7 + 8 ) region towards the 
melting point (Figure 2.1(a)). Both Ae 3 and Ae 1 temperatures are lowered. 

• Type 2 Expanded 7 field - This is same as Type 1 except that sometimes austenite 
range is shortened by the formation of iron-rich compounds as shown in Figure 2.1(b). 
Alloying elements like C, N, Cu, Zn, and Au, belong to this category. The presence 
of C and N expands the 7 field such that the solid solubility of C and N in austenite 
increases to 2.0 wt.% and 2.8 wt.%, respectively. 

• Type 3 Closed 7 field - The alloying elements of this group restrict the temperature 
range for stable austenite. Thus, a smaller 7 loop and a continuous and wider 8 - and 
a-phase fields are obtained as shown in Figure 2.1(c). Alloying elements belonging to 
this category are Si, Al, Be, and the carbide forming elements like Cr, Ti, V, Mo, and 
W. 

• Type 4 Contracted 7 field - In this category, the a- and 7-phase fields are bounded by 
a miscibility gap, i.e., a and 7 solid solutions axe in equilibrium with an inter-metallic 
compound or a solid solution (Figure 2.1(d)). Boron is an important element of this 
group and other elements are the carbide formers like Ta, Zr, and Nb [2]. 

The presence of alloying elements in iron affects the phase diagrams in various ways. In 
order to find their effect on transformation temperature binary iron-alloys phase diagrams 
were studied. Relevant portions of binary phase diagrams of iron, i.e., Fe-C, Fe-Mn, Fe-Si, 
Fe-Ni, Fe-Cr, Fe-Mo, and Fe-Cu are given in Figures 2.2 to 2.8, respectively. The phase 
diagrams are taken from Kubachewski [13]. 
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Figure 2.2: Fe-C Phase diagram 
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Figure 2.3: Fe-Mn Phase diagram 
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Figure 2.4: Fe-Si Phase diagram 
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Figure 2.8: Fe-Cu Phase diagram 
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2.2 Effect of Alloying Elements on Fe-FegC Phase Di- 
agrams 

The presence of alloying elements changes the critical range of phase stability, the position of 
eutectoid point, and the location of a- and 7-phase fields as indicated by the binary Fe-Fe3C 
diagram. The austenite stabilizers lower the eutectoid temperature, thus widening the range 
of temperature over which austenite is stable. The ferrite stabilizers raise the eutectoid 
temperature thereby restricting the 7-phase field. Figure 2.9 shows the effect of alloying 
elements on eutectoid temperature. By proper combination of alloying elements, austenite 
phase can be obtained at room temperature; as in austenitic stainless steel. Eutectoid 
carbon composition also changes with presence of alloying elements. All alloying elements 
tend to decrease the eutectoid carbon content (Figure 2.10) [3]. Usually, no new structural 
constituents appear, although the composition and proportion of phases present in the steel 
may be altered considerably [4]. In certain cases, new structural constituents appear, e.g., 
when Cu exceeds 1% by weight, elemental Cu appears as a second phase. 

2.3 Available Data on Ae 3 Temperature 

Experimental data for 109 multicomponent steels were collected from different sources. Ta- 
ble 4.2 gives the compositions of steels and their experimental Ae 3 temperature. The first 
nineteen steels (SI. No. 1-19) were reported by Grange [5]. The steels were isothermally 
transformed at a number of temperatures for varying periods extending upto several hours 
such that the amount of austenite did not change appreciably even by doubling the holding 
time. The temperature at which the last trace of ferrite transformed to austenite was taken 
as the Ae 3 temperature. 

The United States Steels Corporation [6] has published isothermal transformation dia- 
grams for several steels. Next sixtyseven steels (SI. No. 20-86) were taken from this source ; 
experimental technique used for determining At 3 temperature was same as Grange [5]. The 
Ae 3 temperature and chemical compositions are taken from the table of Baganis [7]. 
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Alloying element, wt*7r 

Figure 2.9: Effect of alloying elements on eutectoid temperature 

Aaranson and Domian [8] austenitized steels samples (SI. No. 87-90) at 1300°C for 15 
minutes and quench them to get martensite. The steels were then reheated to temperatures 
in the vicinity of suspected Ae 3 temperatures and soaked for 15 minutes. The starting 
structure was fine unstable martensite. So, 15 minutes holding temperature was considered 
to be adequate for the establishment of equilibrium. The Ae 3 temperature was taken to lie 
between the highest temperature at which ferrite remained in the micro-structure and the 
lowest temperature at which it was not present. 

Hall et al. [9] homogenized steel samples (SI. No. 91-95) at 1300°C for three days. 
The samples from these steels were then austenitized at 1300°C for 15 minutes and held 
isothermally up to 24 hours. 

Gilmour et al. [10] reported the data for steels (SI. No. 96-98). Samples were homogenized 
at 1200°C for 18 hours and then quenched in water to get a fully martensitic structure. The 
samples were then reheated to the holding temperatures, isothermally transformed for three 
days and quenched in water. Austenite in equilibrium with ferrite was considered to have 
the holding temperature as its Aez temperature. Austenite composition was determined by 
micro-probe analysis. 
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Figure 2.10: Effect of alloying elements on eutectoid carbon content 

Kirchner and Uhrenius [11] allowed steel samples (SI. No. 99-105), which were primarily 
Fe-Cr-Mn, to isothermally transform during one week to seven months. Austenite in equi- 
librium with ferrite at a particular isothermal transformation temperature was considered 
to have the isothermal temperature as its Ae 3 temperature. The composition in this case 
also was determined by. micro-probe analysis. Carbon content of these steels was less than 
0.004% and hence assumed to be zero. 

Swinden and Woodhead[12] austenitized four Fe-C-Cr steel samples (SI. No. 106 109) 
at 1000°C and then allowed isothermal transformation to take place for several hours. 

The Ae 3 temperature for binary alloys were taken from the binary phase diagrams given 

in Kubachewski [13]. 

Most of the data sources did not report the experimental and instrumental errors sep- 
arately. Grange [5] had rounded off the results by a maximum of 5°F i.e. 2.8°C. So the 
error level can be considered to be lower than 2.8°C. No error values were reported in the 
isothermal transformation diagrams [6]. Aaranson and Domian [8] reported an error value 
of approximately ±3°C in their data. Hall et al. [9] reported the experimental error to be 
within ±2°C and Gilmour et al. [10] reported it to be approximately ±1°C. Kirchner et 




CHAPTER 2. LITERATURE REVIEW 


14 


al. [11] reported an error less than ±0.4°C. 

2.4 Existing Models for Predicting Ae 3 Temperature 

Grange [5] analysed the experimental data for 19 steels containing C, Mn, Si, Ni, Cr, and 
Mo and suggested the following model for prediction of Ae 3 temperature : 

T e (in°C) = 854.4 - 179AX C - 13.9X Mn + UAX Si - 17.8X m - 1.7X Cr (2.1) 


where A" represents the weight % of the alloying elements. 

Andrews [14] improved the above model by considering few more data and incorporating 
more alloying elements like Cu, N, V, and P. The equation is : 

T e (in°C) = 9l3-AT-25X Mn -nXcr-20X Cu + 60Xsi + 60X Mo 
+40A jv T lOOAy + 700Ap 

where AT accounts for the effects of C and N on Ae 3 temperature. AT depends on equivalent 
carbon content which is obtained by adding one tenth of the weight % of Ni to the weight 
% of C. The values of AT for different equivalent C contents are given in Table 2.1. 

Table 2.1: Values of AT for Various Equivalent Carbon Content 


equivalent 

C content 

0.05 

0.1 

0.15 

0.2 

0.25 

0.3 

0.35 

0.4 

0.45 

0.5 

0.6 

0.7 

AT(°C) 

24 

48 

64 

80 

93 

106 

117 

128 

137 

145 




Baganis [7] calculated the Ae 3 temperature for different steels using the published data 
on free energy and activity of binary and ternary alloys of iron. According to him, the 
equation is : 

T e = T 0 + AT (2.3) 


where, 


AT = RT 2 


i=2 

[C x D] 
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exp 


' AGijTp) 

RT 0 


+ EuXi n 


B — (1 + euXi tl )exp 


'AGg(To) e u 


RT 0 2 
C = (1 - X li7 ) A H° 0 


(*i, 7 ) 5 


D — exp 


( AGg(r 0 ) e lw 2 ' 
V i?T 0 2 ^ 1 ’ 7 ' , 


T 0 

X li7 

AG? (T 0 ) 

Eu 

AH° 0 


Ae 3 temperature of plain carbon steel, 
mole fraction of the i th alloying element in 7. 

change in the partial molar Gibbs free energy of the i th transformation 
element at temperature To from ferrite to austenite, 
interaction coefficient for the solute i and carbon in iron, 
change in the partial molar enthalpy of iron for ferrite to austenite 
transformation. 

In the expression for T e , the suffixes 0, 1, . . . , i , . . . , n of G, H , and e indicate the con- 
stituents, O stands for iron, 1 for carbon and 2 to n for other substitutional alloying elements. 

Equations 2.1 and 2.2 contain only linear terms which were obtained by fitting straight 
lines to the individual binary systems. Khaira [15], Manoj [16] and Arnab [17] considered 
the self- and mutual-interaction terms. This resulted in considerable decrease in the error in 
the prediction of Ae 3 temperature. According to Khaira [15] : 


T e = 910 - 326.73X C - 32.68X Mn + 45.21X 5i - 36.26*™ + 44.11X c X M n 

+69.74X C X™ - 6.82*£ r + 23.13*^ _ 74.83X c Xf, + 95.69X1, 4) 

-9.46X C X£, + 214.43X£ - 46.05X1, + O.llX^, + 0.17X£ r 
-74.46X C X^ 0 + 37.37X2 X2 r + 8A8 X 2 Ni X C rX Mo 

According to Manoj [16] : 

T e = 912-423.0Xc-36.83X A , n + 93.08X5i-35.68X™~16.78Xcr + 22.17X A * o 

-23.618X Cu + 247.6X£ + 83.56X c X Mn - 124.3X c X Si + 52.44X C X™ 

+41.55XqXcr — 64.09X^Xj^ n 4- 8.541X^ n X™ 4- 436.4Xs,Xc u 

+5.503X^ o — 13.36XcXM n X™ 4- 10.05XcXjtf n X™X<> — 64.66 XcXj^ 0 

+0.8889XAf n X™ — lU.lX Mn X S iX m X Mo + 0.0992Xjw n X^- r -I- 8.971 X^XcyX^o 

(2.5) 
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According to Arnab [17] : 

T e = 912 - 660 X° 88 + 454 11 - 38X™ 2 + 10A|f + 71X\? - 516A^ 18 
+468 A^ 22 - 196 Ac' 84 + 180A^f + 1.13A 4 *- 99 + 29.4A^ 24 - 24A 
+38. 07 Ac 4 A^ 4 + 58. 48 Ac' 33 A^ 46 - 78.89A° 32 Ajf 4 - 6.54A£ 33 Aj^ 7 
-68.9A $ 47 Afc*> + 39. 64 Ac 30 Ac' 86 + 4.48A^ 88 A^ 31 A^ 0 9 

The terms A£, Af^, etc., represent self-interaction (viz., C-C interaction and Mo-Mo 
interaction) and terms like XcXm 0 represent mutual-interaction between the elements, where 
X represents the weight % of the alloying elements. Equations 2.4 and 2.5 were obtained by 
using stepwise multiple linear regression programs given in IMSL. Equation 2.6 was obtained 
by using genetic algorithm. 

The analysis of results obtained from Equations 2.1 to 2.6 on more than 160 steels of dif- 
ferent compositions shows that the standard error is the largest(33°C) in case of Grange’s [5] 
model (Equation 2.1) and the least(8.9°C) in case of Manoj’s [16] model (Equation 2.5). 
The coefficient of correlation for Equation 2.6 is approximately 0.991. Equations 2.4 and 2.5 
do not converge to the respective binary systems, i.e., when all alloying elements except C 
are zero, the equation does not reduce to Ae 3 temperature equation of binary Fe-C system. 
Arnab [17] avoided this problem by forcing binary coefficients into multicomponent model 
(Equation 2.6). 

2.5 Curve Fitting Technique 

Curve fitting technique is used to find the influence of independent variables on dependent 
variable. Regression analysis is the most common tool to handle such kind of problems. 
The relationship between two variables can be represented by simple algebraic expressions, 
e.g., polynomials, exponential or logarithmic functions, known as regression function. The 
regression function involves a set of unknown parameters. Values of parameters, which 
give the best fit for a given set of data can be obtained by regression analysis. Thus, the 
relationship between dependent and independent variable is obtained [18, 19]. 
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2.5.1 Linear Regression 

When the regression model is linear in parameters (but not necessarily in the independent 
variables), then the regression model is a linear regression model. Otherwise the model will 
be a nonlinear regression model [19]. 

Let a dependent variable y be influenced by p independent variables Xi,X 2 , ■ ■ . ,x p and 
has p 4- 1 parameters b 0 , b u . . . , b p . Then, the functional relation is : 

y f (b 0 , bi , . . . , b m , Xi , ^2) • • ■ , -£p) (2.7) 

There are three possibilities : 

1. y is linear with respect to both independent variables and parameters, e.g., 

y — bo + b\X\ 4- ^22-2 + • • • + bpXp (2-8) 

2. y is linear with respect to parameters, but nonlinear with respect to the independent 
variables, e.g., 

y — bo + b\x + bix* + • • • 4- b p x p (2-9) 

3. y is nonlinear with respect to both parameters and independent variables, e.g., 

y = ax b + cx d (2-10) 

The first two equations represent linear regression model where as the third equation 
belongs to nonlinear category. A simple linear regression model is : 

Vi = 0o + fax* + i = l,2,...,n (2.11) 

where the observed values of yi 1 s constitute the dependent variable, the Xi% are the inde- 
pendent variables, /3o and Pi are the intercept and slope parameters, respectively, and e^s 
are independently distributed normal errors, each with mean zero and variance a 2 . 

2.5.2 Multiple Linear Regression 

In multiple linear regression, the dependent variable is regressed on k(k > 1) independent 
variables. A multiple regression model is : 

yi = P 0 + Pl X i\ + @2 x i2 H h Pk x ik + £i t = 1, 2, . . . , H (2.12) 
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where yP s represents the dependent variable, x^’s, x^'s, . . . , x^'s are the settings of the k 
independent variables, /?o, Pi, ■ ■ ■ ,Pk are the regression coefficients and e t -’s are independently 
distributed normal errors with mean zero and variance a 2 [18,19,21]. 

After deciding the mathematical model, the parameters in the model have to be esti- 
mated. Several criteria can be used to estimate the parameters. Least square method is 
frequently used for this purpose. In least square method, the parameters for which the sum 
of the squares of the errors (SSE) is minimum, are obtained. 

SSE = '£J = ’L(Y i -Y i ) 2 (2.13) 

1=1 i= 1 

where, YJ = observed value, and 

V = predicted value of Y{ by the equation. 

2.5.3 Stepwise Regression 

Let Equation 2.12 be the regression model, and z\ , z^, . . . , z T all be function of one or more of 
the x’s, representing the complete set of variables from which the equation is to be choosen. 
Two opposite criteria for selecting a resultant equation are usually involved : 

1. To make the equation useful for predictive purposes the model should include as many 
z’s as possible, so that reliable fitted values can be obtained. 

2. To make the evaluation easy, the equation should include as few z’s as possible. 

The stepwise selection procedure involves insertion of variables until the regression equa- 
tion is satisfactory. The order of insertion is determined by using the partial correlation 
coefficient as a measure of the importance of variables not yet in the equation. Initially, the 
z having the highest correlation coefficient with Y is entered into the equation. The signifi- 
cance of the variable is then tested. In the next step, the partial correlation coefficients of all 
the predictors not in the models is calculated, and the z having the highest value is entered 
into the model. The overall regression is now checked for the significance, the improvement 
in R 2 value is noted and the partial F-values for both the variables in the equation are ex- 
amined. The lower of the two partial F’s is then compared, and the corresponding predictor 
variable is retained in the equation or rejected, according to whether the test is significant 
or not. A predictor that may have been the best entry candidate at an earlier stage, may 
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at stage, be superfluous because of the relationships between it and other variables, now in 
the regression. To check on this, the partial F criterion for each variable in the regression 
at any stage of calculation is evaluated. The process of adding and dropping of predictors 
continues till no further addition or dropping of variables into or from the model, occurs. 

The q; values used for the partial F-test, represents the significance level. One can select 
different values of a for entering into or leaving from the model. a en t ry is always kept less 
than a exit , and it is usually the best to set a entry = a exit [18]. 


2.6 Scope of Present Work 

All models discussed in section 2.4 have certain drawbacks. This could be due to certain 
assumptions in deriving the models. For example, the linear dependence of Ae 3 tempera- 
ture on composition by Grange [5] ignored both self- and mutual-interaction of the alloying 
elements. Andrews [14] improved the Grange’s model [5] by considering self-interaction of C 
and Ni and mutual-interaction between them. But, the interaction terms of other alloying 
elements were ignored. His model was not accurate for steels having Mn>1.0% and Si>0.8% 
by weight, as the interaction terms become important in these steels. Baganis [7] considered 
the equilibrium between ferrite and austenite in his model and used thermodynamic data of 
binary and ternary alloys of iron for calculation. Self-interaction of C and its interaction with 
other alloying elements were taken into account in his model. Self- and mutual-interactions 
of alloying elements other than C, were ignored. His model failed to predict the Ae 3 tem- 
perature, in case of high alloy steels. Baganis’s model (Equation 2.3) also suffers from the 
criticism that it is very complicated and that it involves too many thermodynamic data 
which are not readily available. 

Some of the above drawbacks were overcome in Khaira’s [15] and Manoj’s [16] models 
(Equations 2.4 and 2.5). They considered all possible self- and mutual-interactions among 
the alloying elements. In their model, the root mean square error was reduced to 8.9°C, which 
was considerably low in comparison to Grange’s [5] and Andrew’s [14] models (Equations 2.1 
and 2.2). But, Khaira’s [15] and Manoj’s [16] model do not predict the Ae 3 temperature of 
the binary iron alloys. On the other hand, Arnab’s [17] model (Equation 2.6) predicts the Ae 3 
temperature quite well for the same system, but only within a small range of composition. 
The root mean square error of 9.5°C as associated with his model is marginally higher than 
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that of Manoj’s [16] and Khaira’s [15] models. 

The present work is carried out with a view to develop a model that would predict the 
Ae 3 temperature -very close to the experimental value. In this context, a polynomial model 
for the multicomponent system is assumed. The effect of seven alloying elements, namely 
C, Mn, Si, Ni, Cr, Mo, and Cu, on Ae 3 temperature of steel is studied individually. Seven 
binary models are obtained considering the composition of each alloying elements over a wide 
range. These binary models are then incorporated in the polynomial model, replacing the 
binary terms present in it. The constant term takes the value of transformation temperature 
for pure iron. The remaining mutually interacting terms in the multicomponent model are 
determined by stepwise multiple linear regression analysis. The performance of this model 
is found to be statistically better than that of existing models in the literature. 



Chapter 3 


MODEL FOR BINARY ALLOYS OF 
IRON 

3.1 Selection of Binary Data and Regression Model 

The binary data is taken from the binary phase diagrams given in Kubachewski [13]. 
Manoj [16] and Arnab [17] calculated the binary iron-alloys equations considering alloying 
elements over a small range of compositions. The nonlinear models developed by Arnab [17] 
have two advantages over linear models developed by Manoj [16], namely, greater accuracy 
in calculation and lesser number of terms in the model. But, these models do not give satis- 
factory result when the composition of an alloying element goes beyond the range over which 
the equation was developed. Since, the binary models are forced into the multicomponent 
model, they should give satisfactory result over the entire range of composition for a particu- 
lar alloying element present in multicomponent steels. Hence, in the present work the binary 
equations are recalculated by taking the binary data over the entire range of composition, 
present in multicomponent steels. Table 3.1 shows the maximum and minimum composition 
of an alloying element in multicomponent steels and binary alloys. 

The nonlinear model used for all the binary systems is : 

Y = o 0 + flix 61 + a 2 x h (3.1) 

where, Y = dependent variable 
x = independent variable 

and the values of a 0 , a u a 2 , h, and b 2 can be determined by doing nonlinear regression on 
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Table 3.1: Range of compositions of various alloying elements 


Elements 

C 

Mn 

Si 

Ni 

Cr 

Mo 

Cu 

multicomponent 

min. comp.(in wt.%) 

0 

0 

0 

0 

0 

0 

0 

steels 

max. comp. (in wt.%) 

0.68 

5.01 

1.83 

5.00 

5.70 

0.83 

1.49 

binary 

min. comp. (in wt.%) 

0 

0 

0 

0 

0 

0 

0 

steels 

max. comp. (in wt.%) 

0.77 

5.91 

1.53 

5.77 

6.54 

2.72 

3.40 


known values of x and Y data [22]. Various statistical analysis are then carried out for the 
evaluation of results. Some of the statistical parameters used in this contexts are described 
below : 


• Standard error of estimate : It is also called as standard deviation of estimate and 
is defined as - 

"t(Vi - Yi) 2 ' 


SE = 


n 


(3.2) 


• Residue of error estimate : It is the difference in the experimental and predicted 
values of a dependent variable. 


Residue —Yi — Yi (3.3) 


• Sum of the squared residue : It is the sum of the square of the errors or residues. 

SSE = f^{Yi - Yf (3.4) 


3.2 Determination of Analytical Expression for Ae 3 Tem- 
perature in Binary Iron- Alloy Systems 

The equation of Ae 3 temperature for binary iron-alloy system is represented as : 

T e = a 0 + ai£ bl + a 2 x h (3.5) 

where, T e — experimental Ae 3 temperature 

x = weight fraction of alloying element 
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In each model, a<j is taken to be 912°C, the austenite to ferrite transformation temperature 
and, the values of ax, a 2 , &i, and 62 are determined by using a standard routine (E04FCF) 
available in NAG FORTRAN library. The above model involves four unknowns, therefore, 
instead of obtaining a unique solution, a set of solutions are obtained. Each set of the 
solutions is again tested by using a standard curve fitting software package PEAKFIT. The 
result having best fit and minimum standard error is accepted. In some cases, e.g., Fe-Mn, 
Fe-Ni, and Fe-Cu, the values of 02 and 62 are found to be very low, and hence they are 
neglected. By doing so, the accuracy of the result is found to be unaffected. Thus, for some 
cases, the equation becomes : 

T e = a 0 + a^ 1 (3.6) 


Table 3 . 2 : The binary models coefficients when x is in weight fraction 


System 

a 0 

ax 

b 1 

a 2 

h 

Fe-C 

912 



75345.25 

1.11 

Fe-Mn 

912 

-1867.0 

0.83 

0 


Fe-Si 

912 

10 9 

3.92 

8509 

1.05 

Fe-Ni 

912 

-884.88 

0.62 

0 

0 

Fe-Cr 

912 

11730.51 

1.70 



Fe-Mo 

912 

1.6 x 10 8 

3.99 

3818 

1.05 

Fe-Cu 

912 

-1415.32 

0.97 

0 

0 


Table 3.2 shows the values of parameters for binary iron-alloys systems, when x is ex- 
pressed in weight fraction of the alloying elements. Table 3.3 shows these values, when x 
is expressed in weight % of the alloying elements. Table 3.4 shows the experimental and 
predicted Ae 3 temperature along with the composition of alloying elements and residues. 
The residue is found to be lesser than that in the existing binary models. 
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Table 3.3: The binary model coefficients when x is in weight % 


System 

do 

(Zi 

b 1 


b 2 

Fe-C 

912 

-660 

0.88 

454 

1.11 

Fe-Mn 

912 

-40.84 

0.83 

0 

0 

Fe-Si 

912 

14.45 

3.92 

67.6 

1.05 

Fe-Ni 

912 

-50.92 

0.62 

0 

0 

Fe-Cr 

912 

4.67 

1.70 

-21.16 

1.14 

Fe-Mo 

912 

1.67 

3.99 

30.0 

1.05 

Fe~Cu 

912 

-16.25 

0.97 

0 

0 


Table 3.4: Predicted Ae 3 temperature and residues in binary models 


System 

Weight% of 
alloying element 

Experimental Ae$ 

temperature 

Predicted Ae 3 

temperature 

Residue 


0.00 

912.00 

912.00 

0.00 


0.05 

882.00 

881.05 

0.95 


0.11 

856.00 

856.56 

-0.56 


0.18 

834.00 

833.73 

0.27 


0.25 

815.00 

814.58 

0.42 


0.33 

795.00 

795.83 

-0.83 

Fe-C 

0.40 

782.00 

781.50 

0.50 


0.47 

769.00 

768.76 

0.24 


0.55 

756.00 

755.81 

0.19 


0.60 

749.00 

748.48 

0.52 


0.66 

740.00 

. 740.38 

-0.38 


0.70 

736.00 

735.37 

0.63 


0.77 

727.00 

727.28 

-0.28 
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Table3.4 : Continued 


System 

Weight% of 
alloying element 

Experimental Ae 3 

temperature 

Predicted Ae 3 

temperature 

Residue 


0.00 

. 912.00 

912.00 

0.00 


0.49 

890.00 

889.41 

0.59 


1.48 

856.00 

855.45 

0.55 


1.97 

840.00 

840.30 

-0.30 

Fe-Mn 

2.46 

826.00 

825.79 

0.21 


2.95 

811.00 

811.76 

-0.76 


3.44 

798.00 

798.12 

-0.12 


3.94 

785.00 

784.55 

0.45 


4.43 

771.00 

771.52 

-0.52 


4.92 

759.00 

758.74 

0.26 


5.42 

746.00 

745.92 

0.08 


5.91 

734.00 

733.56 

0.44 


0.00 

912.00 

912.00 

0.00 


0.13 

920.00 

919.94 

0.06 


0.25 

927.00 

927.83 

-0.83 


0.38 

936.00 

936.80 

-0.80 


0.51 

946.00 

946.37 

-0.37 

Fe-Si 

0.63 

956.00 

955.98 

0.02 


0.76 

968.00 

967.60 

0.40 


0.89 

980.00 

980.97 

-0.97 


1.02 

996.00 

996.64 

-0.64 


1.14 

1013.00 

1013.72 

-0.72 


1.27 

1035.00 

1035.76 

-0.76 


1.53 

1094.00 

1094.19 

-0.19 
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Table3.4 : Continued 


System 

Weight% of 
alloying element 

Experimental Ae 3 

temperature 

Predicted Ae 3 

temperature 

Residue 


0.00 

912.00 

912.00 

0.00 


0.53 

878.00 

877.65 

0.35 


1.05 

860.00 

859.52 

0.48 


1.58 


844.38 

0.62 


2.10 

831.00 

831.34 

-0.34 

Fe-Ni 

2.62 

819.00 

819.48 

-0.48 


3.15 


808.29 

0.71 


3.67 

798.00 

797.98 

0.02 


4.20 

789.00 

788.03 

0.97 


5.24 

770.00 

769.81 

0.19 


5.77 


761.06 

-0.06 


0.00 

912.00 

912.00 

0.00 


0.37 

906.00 

906.05 

-0.05 


0.93 

896.00 

896.65 

-0.65 


1.86 

883.00 

882.48 

0.52 

Fe-Cr 

2.80 

870.00 

870.45 

-0.45 


3.73 

860.00 

860.88 

-0.88 


4.67 

854.00 

853.53 

0.47 


5.61 

849.00 

848.49 

0.51 


6.55 

846.00 

845.69 

0.31 


7.49 

845.00 

845.10 

-0.10 
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Table3.4 : Continued 


System 

Weight% of 

Experimental Ae 3 

Predicted Ae 3 

Residue 


alloying element 

temperature 

temperature 



0.00 

912.00 

912.00 

0.00 


0.34 

921.00 

921.78 

-0.78 


0.69 

933.00 

932.90 

0.10 


1.03 

945.00 

945.13 

-0.13 

Fe-Mo 

1.37 

961.00 

960.03 

0.97 


1.71 

979.00 

979.42 

-0.42 


2.04 

1005.00 

1004.77 

0.23 


2.38 

1040.00 

1040.43 

-0.43 


2.55 

1063.00 

1062.92 

0.08 


2.72 

1089.00 

1089.14 

-0.14 


0.00 

912.00 

912.00 

0.00 


0.26 

907.00 

907.60 

-0.60 


0.53 

903.00 

903.22 

-0.22 


0.92 

897.00 

897.01 

-0.01 


1.05 

895.00 

894.96 

0.04 

Fe-Cu 

1.31 

891.00 

890.88 

0.12 


1.57 

887.00 

886.83 

0.17 


1.84 

883.00 

882.64 

0.36 


2.62 

871.00 

870.64 

0.36 


2.88 

867.00 

866.66 

0.34 


3.40 

858.00 

858.74 

-0.74 


3.3 Assessment of Models 

Table 3.5 shows the values of different statistical parameters described in section 3.1. The 
values indicate excellent reliability of models. The standard error in each case is found to be 
very small. Mean residue in most of the cases lies within ±0.08. Binary models for Fe-Mn, 
Fe-Ni, and Fe-Cu individually contain only two terms, but the other four models have three 
terms each. Small number of terms makes the equation easy to evaluate. In all the models, 
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the absolute value of residue is less than 1.00. A good agreement between exper im ental value 
and predicted value can be seen in Figures 3.1 to 3.7. 

Table 3.5: Statistical analysis of the results in binary models 


System 

No. of data 

points(n) 

SSE 

Mean of Residue 

Maximum 

Residue (absolute) 

SE 

Correlation 

coefficient 

Fe-C 

13 



0.95 

0.51 

0.9999624 

Fe-Mn 

12 

2.12 

0.07 

0.76 

0.42 

0.9999708 

Fe-Si 

12 

4.11 

-0.40 

0.83 

0.59 

0.9999667 

Fe-Ni 

11 

2.57 

0.22 

0.71 

0.48 

0.9999527 

Fe-Cr 

10 



0.88 

0.48 

0.9998191 

Fe-Mo 

10 

2.02 

-0.05 

0.78 

0.45 

0.9999709 

Fe-Cu 

11 



0.74 

0.35 

0.9997717 



Wt. % of C 


Figure 3.1: Experimental data and predicted Ae 3 temperature profile of Fe-C system 
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Wt. % Of Mn 


Figure 3.2: Experimental data and predicted Ae 3 temperature profile of Fe-Mn system 



Figure 3.3: Experimental data and predicted Ae 3 temperature profile of Fe-Si system 
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Figure 3.4: Experimental data and predicted ,4e3 temperature profile of Fe-Ni system 



Figure 3.5: Experimental data and predicted Ae 3 temperature profile of Fe-Cr system 
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Figure 3.6: Experimental data and predicted *4e 3 temperature profile of Fe-Mo system 



Figure 3.7: Experimental data and predicted Ae 3 temperature profile of Fe-Cu system 





Chapter 4 

MODEL FOR MULTICOMPONENT 
SYSTEM 


4.1 Formulation of Model 


Austenite in low alloy steels is a dilute solution of fee iron and several alloying elements. The 
Ae 3 temperature of a steel depends on the composition of these alloying elements present 
in it. Hence, in principle, Ae 3 temperature can be estimated in terms of weight fractions of 
various alloying elements present in steel. In this analysis, the following alloying elements 
are considered : C, Mn, Si, Ni, Cr, Mo, and Cu. 

Consider a steel containing 1 to n alloying elements and the solvent is o. The Ae 3 
temperature (T e ) can be expressed as a function of its composition. 

T e — f (xj , X2, . . . , X { , . . . , x n ) (LI ) 


where .r* is the weight fraction of i th alloying elements. 

Assuming the relation to be a polynomial, Equation 4.1 can be expressed as : 


£-E-X> 




Ap 


yrrii ym n 

• * • * -^n 


(4.2) 


mi 




m n 


where A represents the coefficients of polynomial and each m stands for all positive integers 
including zero. 

By Taylor’s series expansion around the composition X 0 — > 1, and Ai, . . . , Xj, . . . , X m — > 0 
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the coefficient of each term can be expressed as derivatives. So, A can be expressed as : 


A. 




1 

d r T e 

flniil 

.axp . . . dX™' . . . dX?" _ 


(4.3) 


JXo4 1 


where, r = 7] rn z 
i 


Equation 4.3 represents the r th order interaction coefficient. Thus, the Ae 3 temperature 
of steel can be expressed as : 

T e = Ao + ^2 AiXi + AuXf + - ^2 Y2 AikXiXk + . . . (4.4) 

i i i fc, kj^i 

where i and k stand for alloying elements and 4 0 is the .4e 3 temperature of pure iron. 

The number of terms in the model depends on the desired accuracy in the results. The 
higher the number of terms, better is the predicted Ae 3 temperature. If only first order 
terms are retained, Equation 4.4 reduces to : 


T e — Aq 4- A\X\ + - ■ • + AiXi ■+••• + A n X n 


(4.5) 


and if only upto second order terms are retained, it reduces to : 

T e = 4Lo + 4i W 4- • ■ • + AiXi + • • • + A n X n + AuXi + • • • + AuXiXi + • • • . . 

(4.bj 

+.4 ln X 1 A r n + A 22 X 2 + 4.23 A 2 A 3 4 1" b A nn X^ 

Similar expressions for T e containing polynomial terms upto any desired order can be written. 


4.2 Evaluation of the Coefficients of the Model 

For evaluating the multicomponent model, a polynomial upto fourth order is considered. 
There are 330 terms in the model consisting of one zero order (or constant) term, 7 first 
order terms, 28 second order terms, 84 third order terms, and 210 fourth order terms. 
Out of these 330 terms, 28 terms belong to binary systems and rest 301 terms belong to 
multicomponent systems. 

The zero order term or the constant term is taken to be 912°C, the austenite to ferrite 
transformation temperature of pure iron. The 28 binary terms in this multicomponent 
model are replaced by the binary terms obtained earlier from the binary models. Even by 
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this substitution, the number of independent variables involved in the calculation is still 
fairly large (301 variables). Stepwise linear regression technique discussed in section 2.5 is 
used for evaluating the important parameters and their coefficients. The regression analysis 
is carried out on the multicomponent data. 109 data sets involving multicomponent steels 
are used for parameter estimation. 

The number of variables retained after regression analysis depended on the significant 
levels, namely, aq and »■>, 35 defined in the calculation involved in the stepwise regression 
technique [19, 20, 21]. In the present work, the values of aq and a 2 are taken to be equal in 
accordance with the reasoning given by Draper & Smith [18]. Finally, the coefficients of the 
multicomponent model are evaluated for different sets of aq and a 2 . It is observed that the 
number of variables increased with the values of significant levels, (i.e., an and a 2 )- • 

A standard subroutine “RSTEP” available in IMSLibrary is used for evaluating the co- 
efficients. Appendix-A shows the calling program. 

Prior to this, a slightly different approach is tried in which the constants and binary 
terms in the multicomponent model are replaced by the binary values obtained by Arnab 
[17], and then a stepwise linear regression is carried out as discussed earlier for the evaluation 
of multicomponent terms. This model is given Appendix-B. This model when compared to 
the one discussed earlier is found to be statistically inferior. The statistical analysis is also 
presented in Appendix-B. 

4.2.1 Coefficients of the Model 

The regression analysis is carried out for different values of ax and a 2 . The variables and 
their coefficients obtained by regression analysis, setting aq and <x 2 as 0.01, 0.03, and 0.05, 
are shown in Table 4.1. The number of multicomponent terms in the equations are found 
to be 9, 10, and 12, respectively. It is noticed that the 9 variables which are present in the 
model by setting oq = 0.01 = oc 2 are also present in the other two models. But, they contain 
few other terms in addition to those 9 variables. The equations so obtained from regression 
analysis are : 
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• Model for «i = a 2 = 0.01 (Model 1) : 

T e = 912 - 660X^ 88 + 454^ u - 40.84X^ 8 n 3 + 14.45X|f + 67.6A 4 ; 05 
-50.92A°f + 4.67A£ r 70 - 21.16A ( M 4 + 1.67A 3 , 99 + 30.0A^° o 5 
-16.25A 9 ' 97 + 66.09A C A^ + 40.8lA c A Cr + 284.9A Mn A 5i 
+ 13.48A c Af fn - 35.46A Mn 3A Si - 48.24A c A| fo - 17.39 X 2 c X 2 m 
-123 5X c X Mn X 2 Si - 547.9A Mn A Si AL 

• Model for ai = ao = 0.03 (Model 2) : 

T e = 912 - 660Ag- 88 + 454A£ U - 40.84A 9 / n 3 + 14.45Af; 25 + 67.6Ajf 
-50.92A'vP H- 4.67 XqJ° - 21.16A£ r 14 + 1.67 A 3 , 99 + 30.0A^° o 5 
- 16.25 A 9 f + 69.4lA c A jVi + 38.12A c A Cr + 312.7A M „A S ; 

+ 13.17A c A 2 fn - 4ZMXl n X Si - 48.74X c X 3 Mo - l8.14X 2 c X 2 Nl 
- UGBXcX^Xli - 67.5lX Mn X Si X 2 Mo - 547.9A Mn A 5l A 2 , 0 

• Model for at = «■> = 0.05 (Model 3) : 

T, = 912 - C60A'c 88 + 454 11 - 40.84A^ 8 n 3 + 14.45 Af; 25 + 67.6 A^ 05 
-50.92 AftP + 4.67A'cJ° - 21.16A^ 14 + 1.6 7 A 3 , 99 + 30.0A^° o 5 
- 16.25 A^f + 69.31 X c X Ni + 37.79A c A Cr + 312.8A Wn A St 
f 13. 12A'cA’(/„ — 40.98A^ n Asi — 50.33AcA^ o — 18.22AcA^- i 
+8525A|rA Mo X C u ~ 1328A C A W „A| < - 59.58A Mn A 5l A^ 0 
- 548. 5 A MnXsiXho - 1.05 x 10 4 A Mn A 5i A Mo A Cu 
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Table 4.2 : Continued 

SI. No. Weight% of alloying elements 4 e , tem 


C 

0.56 

0.59 

0.37 

0.10 

0.40 

0.38 

0.42 

0.33 

0.11 

0.65 

0.26 

0.35 

0.48 

0.68 

0.42 

0.40 

0.36 

0.22 

0.40 

0.41 

0.39 

0.50 

0.48 

0.38 

0.11 


Weight% of alloying elements 
Mn Si Ni Cr Mo 

0.26 0.21 1.97 0.04 0.00 

0.25 0.21 3.90 0.04 0.01 

0.68 0.21 3.41 0.05 0.00 

0.52 0.28 5.00 0.07 0.03 
0.52 0.28 5.00 0.07 0.03 

0.37 0.25 0.07 0.57 0.00 

0.68 0.16 0.07 0.93 0.00 

0.45 0.28 0.06 1.97 0.02 

0.38 0.44 0.15 5.46 0.42 

0.68 0.23 0.02 0.02 0.10 

0.87 0.29 0.03 0.06 0.26 

0.80 0.24 0.06 0.06 0.25 

0.94 0.16 0.04 0.03 0.25 

0.87 0.26 0.01 0.03 0.24 

0.20 0.18 0.00 0.00 0.21 

0.42 0.19 0.06 0.13 0.53 

0.17 0.16 0.00 0.00 0.82 

0.79 0.24 0.03 0.08 0.50 

0.57 0.03 3.49 0.00 0.01 

0.60 0.03 3.51 0.00 0.21 

0.56 0.03 3.53 0.00 0.74 

0.23 0.53 0.02 0.05 0.00 

0.57 0.20 0.06 0.02 0.01 

0.72 0.21 1.32 0.49 0.00 

0.45 0.18 3.33 1.52 0.03 

0.53 0.27 0.00 0.90 0.18 

0.77 0.15 0.04 0.98 0.21 

0.57 0.27 1.87 0.45 0.24 

0.78 0.24 1.79 0.80 0.33 

1.90 0.05 0.24 
1.84 0.06 0.23 


Aez temperature in °C 

Residue 

Expt . 

Pred. 

in °C 

732.00 

738.78 

- 6.78 

712.00 

701.09 

10.91 

734.00 

742.58 

- 8.58 

774.00 

771.64 

2.36 

714.00 

715.35 

- 1.35 

791.00 

791.22 

- 0.22 

782.00 

773.87 

8.13 

807.00 

802.85 

4.15 

839.00 

846.68 

- 7.68 

741.00 

744.53 

- 3.53 

824.00 

834.21 

- 10.21 

802.00 

801.66 

0.34 

769.00 

765.17 

3.83 

744.00 

732.88 

11.12 

782.00 

790.83 

- 8.83 

802.00 

781.34 

20.66 

800.00 

800.85 

. - 0.85 

814.00 

824.79 

- 10.79 

724.00 

715.29 

8.71 

718.00 

718.07 

- 0.07 

718.00 

725.20 

- 7.20 

780.00 

780.77 

- 0.77 

764.00 

761.85 

2.15 

754.00 

764.35 

- 10.35 

782.00 

764.43 

17.57 

800.00 

815.43 

- 15.43 

786.00 

7 , 85.96 

0.04 

810.00 

802.67 

7.33 
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SI. No. 

Weight.% o 

— — * t *— — r* — — 


C Mn 

Si 

62 

0.16 0.52 

0.27 

63 

0.18 0.79 

0.21 

64 

0.30 0.80 

0.29 

65 

0.59 0.89 

0.24 

66 

0.44 .90 

0.25 

67 

0.24 0.94 

0.47 

68 

0.10 1.63 

0.21 

69 

0.63 0.87 

0.17 

70 

0.14 0.81 

0.34 

71 

0.16 0.60 

0.27 

72 

0.18 0.57 

0.24 

73 

0.43 1 .02 

0.26 

74 

0.45 0.89 

0.31 

75 

0,50 0.77 

0.23 

76 

0.19 0.77 

0,24 

i < 

0.18 0,67 

: 0,17 

78 

0.19 0.75 

0.18 

79 

0.17 1.65 

0.16 

80 

0.21 0.75 

0.18 

81 

0.21 0.78 

0.19 

82 

0.22 0,77 

0.19 

83 

0.18 0.65 

0.12 

84 

0.18 0.75 

0.71 

85 

0.37 0.77 

0.15 


Ae 3 tempera ture in °C 
Expt. Pred. 


0.15 I 0.04 I 0.98 0.21 
0.08 0.50 
2.99 0.00 
2.99 0.00 
0.00 0.39 
0.00 0.38 
0.00 0.00 
0.00 0.00 


777.00 

838.00 

792.00 

745.00 

775.00 

844.00 

864.00 

736.00 

810.00 

800.00 

850.00 

770.00 

773.00 

755.00 

827.00 

824.00 

823.00 

800.00 

811.00 

800.00 

786.00 

811.00 

855.00 

787.00 


Pred. 

784.54 

821.57 

800.39 

742.28 
767.89 

843.21 
849.99 
742.51 
817.84 
807.60 

835.28 
773.77 
766.03 

756.83 

829.40 

803.22 
803.44 

789.19 

793.84 
791.15 

783.19 
796.74 
853.39 
785.96 


Residue 
in °C 
- 7.54 
16.43 
- 8.39 

2.72 
7.11 
0.79 
14.01 
- 6.51 
- 7.84 
- 7.60 

14.72 
- 3.77 
6.97 
- 1.83 
- 2.40 
20.78 
19.56 
10.81 
17.16 
8.85 
2.81 
14.26 
1.61 
1.04 









mperature in °C Residue 


Pred. 

772.89 
741.10 
902.01 

722.89 

739.55 
740.51 
802.78 

807.00 
789.46 
813.70 
738.93 
748.15 
738.09 

817.72 
811.67 
814.63 

812.00 

912.00 
881.05 

856.56 

833.73 
814.58 
795.83 
781.50 
768.76 
755.81 
748.48 
740.38 
735.37 
727.28 
889.41 


in °C 


21.49 


10.54 

- 13.70 

11.07 


11.91 

11.28 


- 8.00 


- 0.83 
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Table 4.2 : Continued 


SI. No . 


Weig 

; ht % of alloying elements 

Ae 3 temperature in °C 

Residue 

in °C 

C 

Mn 

Si 

Ni 

Cr 

Mo 

Cu 

Expt. 

Pred. 

124 

0.00 

1.48 

0.00 

0.00 

0.00 

0.00 

0.00 

856.00 

855.45 

0.55 

125 

0.00 

1.97 

0.00 

0.00 

0.00 

0.00 

0.00 

840.00 

840.30 

-0.30 

12 G 

0.00 

2.46 

0.00 

0.00 

0.00 

0.00 

0.00 

826.00 

825.79 

0.21 

127 

0.00 

2.95 

0.00 

0.00 

0.00 

0.00 

0.00 

811.00 

811.76 

-0.76 

128 

0.00 

3.44 

0.00 

0.00 

0.00 

0.00 

0.00 

798.00 

798.12 

-0.12 

129 

0.00 

3.94 

0.00 

0.00 

0.00 

0.00 

0.00 

785.00 

784.55 

0.45 

130 

0.00 

4.43 

0.00 

0.00 

0.00 

0.00 

0.00 

771.00 

771.52 

-0.52 

131 

0.00 

4.92 

0.00 

0.00 

0.00 

0.00 

0.00 

759.00 

758.74 

0.26 

132 

0.00 

5.42 

0.00 

0.00 

0.00 

0.00 

0.00 

746.00 

745.92 

0.08 

133 

0.00 

5.91 

0.00 

0.00 

0.00 

0.00 

0.00 

734.00 

733.56 

0.44 

134 

0.00 

0.00 

0.13 

0.00 

0.00 

0.00 

0.00 

920.00 

919.94 

0.06 

135 

0.00 

0.00 

0.25 

0.00 

0.00 

0.00 

0.00 

927.00 

927.83 

-0.83 

13 C 

0.00 

0.00 

0,38 

0.00 

0.00 

0.00 

0.00 

936.00 

936.80 

-0.80 

137 

0.00 

0.00 

0.51 

0.00 

0.00 

0.00 

0.00 

946.00 

946.37 

-0.37 

138 

0.00 

1 

0.00 

0.63 

0.00 

0.00 

0.00 

0.00 

956.00 

955.98 

0.02 

130 

! 0.00 

0.00 

0.76 

0.00 

0.00 

0.00 

0.00 

968.00 

967.60 

0.40 

140 

0.00 

! o.oo 

0.89 

0.00 

0.00 

0.00 

0.00 

980.00 

980.97 

-0.97 

141 

0.00 

0.00 

1.02 

! o.oo 

0.00 

0.00 

0.00 

996.00 

996.64 

-0.64 

142 

0.00 

0.00 

1.14 

0.00 

0.00 

0.00 

0.00 

1013.00 

1013.72 

-0.72 

143 

0.00 

0.00 

1.27 

0.00 

0.00 

0.00 

0.00 

1035.00 

1035.76 

-0.76 

144 

0.00 

0.00 

1.53 

0.00 

0.00 

0.00 

0.00 

1094.00 

1094.19 

-0.19 

145 

0.00 

0.00 

0.00 

0.53 

0.00 

0.00 

0.00 

878.00 

877.65 

0.35 

146 

0.00 

o.oo 

0.00 

1.05 

0.00 

0.00 

0.00 

860.00 

859.52 

0.48 

147 

0.00 

0.00 

0.00 

1.58 

0.00 

0.00 

0.00 

845.00 

844.38 

0.62 

148 

0.00 

0.00 

0.00 

2.10 

0.00 

0.00 

0.00 

831.00 

831.34 

-0.34 

149 

0.00 

0.00 

0.00 

2.62 

0.00 

0.00 

0.00 

819.00 

819.48 

-0.48 

150 

0.00 

0.00 

0.00 

3.15 

0.00 

0.00 

0.00 

809.00 

808.29 

0.71 

151 

0.00 

0.00 

0.00 

3.67 

0.00 

0.00 

0.00 

798.00 

797.98 

0.02 

152 

0.00 

0.00 

0.00 

4.20 

0.00 

0.00 

0.00 

789.00 

788.03 

0.97 

153 

0.00 

0.00 

0.00 

5.24 

0.00 

0.00 

0.00 

770.00 

769.81 

0.19 

154 

0.00 

0.00 

0.00 

5.77 

0.00 

0.00 

0.00 

761.00 

761.06 

-0.06 
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SI. No. Weight % of alloying elements 

C Mn Si | Ni | Cr | Mb" 

155 0.00 0.00 0.00 0.00 0.37 0.00 

156 0.00 0.00 0.00 0.00 0.93 0.00 

157 0.00 0.00 0.00 0.00 1.86 0.00 

158 0.00 0.00 0.00 0.00 2.80 0.00 

159 0.00 0.00 0.00 0.00 3.73 0.00 

160 0.00 0.00 0.00 0.00 4.67 0.00 

161 o.oo 0.00 0.00 0.00 5.61 0.00 

162 0.00 0.00 0.00 0.00 6.55 0.00 

163 0.00 0.00 0.00 0.00 7.49 0.00 

164 0.00 0.00 0.00 0.00 0.00 0.34 

165 0.00 0,00 0,00 0.00 0.00 0.69 


166 0.00 0.00 0.00 0.00 0.00 1.03 

167 0.00 | 0.00 0.00 0.00 0.00 1.37 

168 0.00 j 0.00 0.00 0,00 0.00 1.71 

169 0.00 | 0.00 0.00 0.00 0.00 2.04 

170 0.00 0,00 j 0,00 0.00 0.00 2.38 

171 0.00 0.00 0.00 0.00 0.00 2.55 

172 0,00 0.00 0.00 0.00 0.00 2.72 

173 0.00 0.00 0.00 0.00 0.00 0.00 

174 0.00 0.00 0,00 0.00 0.00 0.00 

175 0,00 0.00 0.00 0.00 0.00 0.00 

176 0.00 0.00 0.00 0.00 0.00 0.00 

177 0.00 0.00 0.00 0.00 0.00 0.00 

178 0.00 0.00 0,00 0.00 0.00 0.00 

179 0.00 0,00 0.00 0.00 0.00 0.00 


Ae^ temperature in ° C Residue 
Expt. Pred. in °C 


906.00 

896.00 

883.00 

870.00 

860.00 

854.00 

849.00 

846.00 

845.00 

921.00 

933.00 

945.00 

961.00 

979.00 

1005.00 


180 0.00 

181 0.00 


0.00 0.00 0.00 0.00 0.00 

0.00 0.00 0.00 0.00 0.00 


182 1 o oo I o no 1 n on I o.oo o.oo I o.oo 


907.00 

903.00 

897.00 

895.00 

891.00 

887.00 

883.00 

871.00 

867.00 

858.00 


906.05 

896.65 

882.48 
870.45 
860.88 
853.53 

848.49 
845.69 
845.10 
921.78 
932.90 
945.13 
960.03 
979.42 
1004.77 


1040.00 1040.43 

1063.00 1062.92 

1089.00 1089.14 


907.60 

903.22 

897.01 

894.96 

890.88 

886.83 

882.64 

870.64 
866.66 
858.74 
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Table 4.3: R 2 and R 2 Adj Values for three models 



R 2 values 

R 2 Adj values 

Model 1 

0.9901 

0.9889 

Model 2 

0.9910 

0.9896 

Model 3 

0.9917 

0.9905 


4.3 Critical Appraisal of Regression Models 

The various models corresponding to different, o, and a 2 values are statistically evaluated 
for the best tit. The results of the statistical analysis for the three models are shown in 
Table 4.4. All the models show comparable statistical parameters. Hence, to select the best 
model out of these three models, additional statistical analysis, R 2 and R 2 Adj and F-test, are 

carried out. 


4.3.1 Comparison of f ? 2 and R~ uij Values 

The values of If and R% () are given by the following expressions. 

H “ 1 " TSS ~ ZW - Y) 2 


(4.10) 


1 


£<;/(N -A') 


where 


RSS 

TSS 

N 

K 
Y 


v Mj “ 1 " _ Y) 2 j{N - 1) 

Residual sum of square 


(4.11) 


True sum of square 
Number of data points 
Number of parameters 

From Table 4.3 and Table 4,4, it is observed that model 1 is statistically inferior to other 
two models, and hence rejected. The R % and ^ues for model 2 and model 3 are nearly 
equal. So, to find the better one among these two models, F-test is conducted. 
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4.3.2 F-test for finding the best model 

The F distribution is a continuous probability distribution. It has two parameters, v x and i/ 2 , 
which are positive whole numbers. The F distribution of two independent variables having 
i y i and U '> . as the degrees of freedom, respectively, is given by : 

(412) 

V\ v 2 

where X{vi) 2 and X{u 2 ) 2 are independent variables. 

The values of F distribution have been tabulated in standard statistical books, for dif- 
ferent values of v\ and v 2 at certain degrees of accuracy, a (also known as confidence level). 
Again 

F(i a) — ~p 

r {v\ ,V2,a) 

F-test can be done in term of number of restrictions. On simplification of the F* expres- 
sion the equation becomes : 

F* — ~ 

1 -RI R /{N-K) 

where m = Number of restrictions 

K = Number of parameters in unrestricted model Comparing the model 2 and 
N — Number of data points 

model 3, it is observed that model 3 contain two extra terms than model 2. So, eliminating 
the extra terms in model 3, model 2 can be obtained.So, model 3 is the unrestricted equation 
and model 2 is the restricted equation, and the number of restrictions is two. To carry out 
the F-test for the best models, out of these two, following assumptions are made : 

• H 0 : Both the models, i.e., model 2 and model 3, are statistically same 
(null hypothesis) 

• H a ■ Model 3 is statistically better than model 2 
(alternate hypothesis) 

And if F(rn,N-K,a) > F*; conclude H 0 
< F *; conclude H a 

Assuming 0.95 confidence level, the value of F* is found to be 8.55, and the value of F( m ^-K,a ) 
i.e., F( 2 ,i 59 ; o. 95 ) from the F distribution table is found to be 19.5, which is greater than 8.55. 
This shows that both the models perform well. 
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Model 3 contains two extra terms than that of model 2. Both the extra terms are of 
fourth order and contain Xq u ■ The fourth order interaction coefficients for Cu are negligible. 
Model 2 containing lesser number of terms, makes the evaluation of Ae :i temperature com- 
putationally easier, hence, accepted as the best model. A better symmetrical distribution of 
residues is also observed in case of model 2. Also the mean of the residues is almost zero for 
model 2. All these reasons goes in favour of selecting the model 2 as the best model. 

Experimental versus predicted Ae 3 temperature curves of steels for the three models are 
shown in Figures 4.1, 4.4, and 4.5. The frequency distribution of residues and cumulative 
frequency distribution of residues are shown in Figures 4.2, 4.3, 4.6, and 4.7. More than 
50% of the residues lie within ±1. Nearly 84% of the residues lie within ±10. All these 
observations enhance the acceptability of the models. 

The error in the model prediction is presumed to be due to the inherent error in the 
experimental data itself. The experimental uncertainties in the measured Ae 3 temperature 
is discussed in section 2.3. Re-examination of the data for which the residues are more than 
±10 is suggested. Chemical analysis should be carried out and the effect of impurities on 
Ae 3 temperature should be critically evaluated. Impurities like phosphorus and sulphur are 
proved to be notorious in this regard. According to Andrews [14], phosphorus influences the 
Ae 3 temperature to the extent of 700°C per weight % of phosphorus. 

4.4 Comparison with Existing Models 

Table 4.3 shows various statistical values of the earlier models and the present model. From 
the table, it can be inferred that standard error in the present model is less compared to the 
existing ones. The correlation coefficient of the present model is also better. The number of 
terms in the present model is 22 which is higher than Arnab’s [17] and Khaira’s [15] models. 
A more symmetric frequency distribution of residue is obtained in the present model. The 
correlation coefficients for all binary systems is more than 0.999 (Table 3.5) for the present 
model. When multicomponent steels are considered, the correlation coefficient is found to be 
0.97 as seen in Table' 4.4. This shows that the model fits both multicomponent and binary 
systems satisfactorily. 



CHAPTER 4. MODEL FOR MULTICOMPONENT SYSTEM 


46 



Experimental Ae 3 temperature (°C ) 


Figure 4.1: Scattergram for model 2 



Cumulative Frequency _ % Frequency Distribution 
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Figure 4.2: Frequency distribution of residues for model 2 



Figure 4.3: Cumulative frequency distribution of residues for model 2 




Predicted Ae 3 temperature (°C ) 
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Figure 4.4: Scattergram for model 1 



Predicted Ae 3 temperature (°C ) 
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Figure 4.5: Scattergram for model 3 
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Figure 4.6: Frequency distribution of residues for model 1 and model 3 



Figure 4.7: Cumulative frequency distribution of residues for model 1 and model 3 
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Table 4.4: Statistical analysis of the results obtained by three models considering both binary 
and multicomponent steels 


Statistical 

Multicomponent steels model for 

parameters 

Oil = 0.01 = Ct 2 

a i = 0.03 = 

c*i = 0.05 = oc2 

No. of steels 

182 

182 

182 

SSE 

10606.13 

9802.77 

8837.09 

SE 

7.64 

7.34 

6.97 

Mean Error 

0.24 

0.0012 

-0.052 

Min. Residue 

-20.78 

-19.43 

-19.83 

Max. Residue 

22.63 

21.49 

21.58 

Correlation 




Coefficients 

0.9950453 

0.9954084 

0.9958596 


Table 4.5: Statistical analysis of the results from different models for .4e 3 temperature 


Statistical 

parameters 

Present 

(model 2) 

Arnab 

[17] 

Manoj 

[16] 

Khaira 

[15] 

Baganis 

[7] 

Andrews 

[14] 

Grange 

No. of steels 

182 

177 



171 

160 

177 

No. of terms 

22 

20 

24 


— 

id 

6 

SE 

7.34 

9.5 

8.9 

9.0 

13.8 

19.2 

33.3 

Residue (max.) 

21.49 

24.88 

34.79 

20.5 

50.1 

85.1 

203.6 

Residue (min.) 

-19.43 

-23.54 

-25.70 

-22.5 

-16.8 

-76.8 

-52.4 

Correlation 

coefficients 


0.9911 




0.9610 

0.909 


SENTBM-Jgr 
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• Fe-Mo : T e = 912 + 1.67X 3 , 99 + 30.0X^° o 5 

. Fe-Cu : T e = 912 - 16.25 X£f 

These correlations are incorporated in a multicomponent model that correlates the At 3 
temperature to the composition of alloying elements. This model unlike the earlier binary 
model contains the mutual interaction terms. A regression analysis of the data led to basically 
three models, each one differing by the levels of significance associated with them. In other 
words, the regression analysis is carried out several times with varying degrees of significance 
and only three models are retained solely on the consideration of number of terms involved in 
the model. The connotation of the significance level is merely a way of keeping and rejecting 
a given variable from the model. 

The model equations are given below: 

• Model for a.\ — c 12 = 0.01 (Model 1) : 

T e = 912 - 660X£ 88 + 454 Ac 11 - 40.84X^ 8 n 3 + 14.45Xff 5 + 67.6X£; 05 
—50.92 Xftf + 4.67XqJ° - 21.16X£ r 14 + 1.67X|/ 0 9 + 30.0X^° 5 
— 16.25-X^' 97 + 66.09XcXAri + 40.81 XqXct + 284.9X;v/ n Xsi 
+13.48X c X& n - 35.46X Mn 3X Si - 48.24X c X 3 Mo - 17.39X 2 X 2 t 
-l235X c X Mn X 2 Si - 547.9 X Mn X s ■>Xlo 

• Model for cti = a 2 = 0.03 (Model 2) : 

T e = 912 - 660X£ 88 + 454X«h 11 - 40.84X^ 83 + 14.45Xf; 25 + 67.6A^ 05 
— 50.92 + 4.67XqJ° - 21.16X£ r 14 + 1.67X|£ 9 + 30.0 A^ 0 o 5 
-16.25X£ 97 + 69.4lA c Ajvi + 38.12X c X Cr + 312.7X Mn X 5i 
+13.17X c X 2 Mn - 43.03X 3 , n X 5i - 48.74X c X 3 Mo - 18.14X 2 c X 2 m 
-1269A c A Mn Al' - 67.5lX Mn X 5i X 2 , 0 - 547.9 X Mn X Si X 2 Mo 

• Model for a\ = c *2 = 0.05 (Model 3) : 

T e = 912 - 660X£ 88 + 454Xq U - 40.84A^ 83 + 14.45X|f + 67.6X^ 05 
-50.92X&P + 4.67A £ r 70 - 21.16X^ r 14 + 1.6 7 A 3 , 99 + 30.0 A^° 0 5 
-16. 25 A 9 f + 69.31 X c X m + 37.79X c X C r + 312.8A M „A Si 
+13.12X c X 2 Mn - 40.98A|, n A 5! - 50.33X c X* Mo - 18.22X 2 c X 2 Ni 
+8525 A^ A MoXcu ~ \328X c X M nX 2 Si - m.38X Mn X Si X 2 Mo 
-54S.5X Mn X Si X 2 Mo - 1.05 x 10 4 X Mn X s *X. Wo X Cu 
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These three models chosen predict the Ae 3 temperatures quite well, both for the binary 
as well as for the multicomponent steels. However, a careful analysis of the above three 
models reveals that in a statistical sense the second model is the most acceptable one. 
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Appendix A 


c ** *** STEPWISE LINEAR REGRESSION ANALYSIS OF *** 

C** *** MULTICOMPONENT STEELS *** 

external rstep,mxtxf,amach 

real d,amach,aov,coef,cov,covs, hist, pin, pout, scale, tol 
integer invoke, iprint,nca,istep, level, nforce, nobs, nstep 

dimension d(109,8),c(109,302),cov(302,302),level(302),scale(302) 
dimension hist(302),aov(13),coef(302,5),covs(302,302) 

nvar=302 

ndata=109 

open(unit = 1, file = ’p.data’) 
read(l,*)((d(n,j),j = 1 , 8), n = 1 , 109) 
do n = 1 , 109 

c(n,302) = d(n,8) - (912.0 - 37979.0358*d(n,l)**0.88 
1 + 75345.245*d(n,l)**l.ll - 1867.0*d(n,2)**0.83 

1 + i e +9*d(n,3)**3.92 + 8509.0*d(n,3)**1.05 
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1 - 884.89*d(n,4)**0.62 + 11730.51*d(n,5)**1.70 

1 - 4031.96*d(n,5)**1.14 + 1.6e+8*d(n,6)**3.99 

1 + 3818.0*d(n,6)**1.05 - 1415.32*d(n,7)**0.97) 

enddo 


C** generating second degree polynomials 

do n = 1 , 109 
m = 1 
do i = 1 , 7 
do j = i+1 , 7 

c(n,m) = d(n,i) * d(n,j) 
m = m + 1 
enddo 
enddo 
enddo 

C** generating third degree polynomials 


20 


£** 


do n = 1 , 109 
m = 22 
do i = 1 , 7 
do 20 j = 1 , 7 

if(j.eq.i) go to 20 
c(n,m) = d(n,j) * d(n,i)**2 
m = m + 1 
continue 
enddo 
enddo 


do n = 1 , 109 
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m = 64 
do i = 1 , 7 
do j = i+1 , 7 
do k = j+1 , 7 

c(n,m) = d(n,k) * d(n,j) * d(n,i) 
m = m+1 
enddo 
enddo 
enddo 
enddo 

C** generating fourth degree polynomials 

do n = 1 , 109 
m = 99 
do i = 1,7 
do 10 j = 1,7 

if(j.eq.i) go to 10 
c(n,m) = d(n,j) * d(n,i)**3 
m = 1 + m 
continue 
enddo 
enddo 


10 




do n = 1 , 109 
m = 141 
do i = 1 , 7 
do j = i+1 , 7 

c(n,m) = d(n,j)**2 * d(n,i)**2 
m = m + 1 


enddo 
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enddo 

enddo 

Q** 

do n = 1 , 109 
m = 162 
do i = 1 , 7 
do 33 j = 1 , 7 
if(j.eq.i) go to 33 
do 34 k = j+1 , 7 
if(k.eq.i) go to 34 

c(n,m) = d(n,k) * d(n,j) * d(n,i)**2 
m = m + 1 

34 continue 

33 continue 

enddo 
enddo 

£** 

do n = 1 , 109 
m = 267 
do i = 1 , 7 
do j = i+1 , 7 
do k = j+1 , 7 
do 1 = k+1 , 7 

c(n,m) = d(n,l)*d(n,k)*d(n,j)*d(n,i) 
m = m+1 
enddo 
enddo 
enddo 
enddo 


enddo 
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read(*,*)a,b 

nca = 302 
nra = 109 

call mxtxf(nra,nca,c,nra,nca,cov,nca) 

invoke = 0 
datalevel /301*1 -1/ 
nforce = 0 
nstep = -1 
istep = 1 
nobs = 109 
pin = a 
pout = b 

tol = 100.0*amach(4) 
iprint = 1 

call RSTEP(invoke, nca, cov, nca, level, nforce, nstep, istep, nobs, pin, 
1 pout, tol, iprint, scale, hist, iend,aov,coef, nca, covs, nca) 

do i = 1,13 

write(*,*)aov(i) 

enddo 

stop 

end 



Appendix B 


During the early stages of model development, Ae 3 temperature of steels are predicted by 
considering a polynomial model as described in the section 4.1 for multicomponent steels. 
For the binary terms, several models available in the literature are compared and Arnab’s [17] 
models are found to be the best amongst all. Accordingly, these binary models are incor- 
porated in the multicomponent model. Fixing the constant terms as 912°C, the important 
multicomponent terms and their coefficients are evaluated using the stepwise linear regression 
over multicomponent data. The equation so obtained is given below : 

T e = 912 - 660X° C 88 + 454A& 11 - 38A^ 9 n 2 + 71A|; 13 - 516Aj^ 8 + 468AM 2 

-19 6A£ 34 + lSOA^ 37 + 1.13Afc°? + 29 AX l { 2 * - 24A°f + 83.95A c X Mn 
+79A8X c X m + 37.72X c X Cr + 305AX Mn X Si - 82.15A£Aa fn (B.l) 

-94.71A& n Asi - 8 93.5Xl.Xc - 13.0 6X 2 Ni X c - 82A5X 2 Mo X c 
-284.3 X c X Mo X Si + 200.1 A 3 , A c + 1029A^A| i - 589.9A| I A Mn A Mo 

Statistical analysis of the results obtained through this equation is performed. The 
maximum and minimum residues are found to be 26.5 and -21.2, respectively. Considering 
all the binary and multicomponent data, the standard error was 8.32. On this basis, the 
model turned out to be statistically better than the existing models. But, it suffers from the 
drawback that it does not predict the Ae 3 temperature satisfactorily in case of binary alloys 
having higher composition. So, in the next attempt the binary equations are developed 
considering the alloying elements over a wide range of composition. 


62 




fYifnS 




a 


m 


l<n r<.-~' nc ’ v 



A125413 



